Wetting and particle adsorption in nanofiows 
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Abstract 

Molecular dynamics simulations are used to study the behavior of closely-fitting spherical 
and ellipsoidal particles moving through a fluid-filled cylinder at nanometer scales. The 
particle, the cylinder wall and the fluid solvent are all treated as atomic systems, and special 
attention is given to the efl'ects of varying the wetting properties of the fluid. Although the 
modification of the solid-fluid interaction leads to significant changes in the microstructure 
of the fiuid, its transport properties are found to be the same as in bulk. Independently of 
the shape and relative size of the particle, we find two distinct regimes as a function of the 
degree of wetting, with a sharp transition between them. In the case of a highly-wetting 
suspending fluid, the particle moves through the cylinder with an average axial velocity 
in agreement with that obtained from the solution of the continuum Stokes equations. In 
contrast, in the case of less-wetting fluids, only the early-time motion of the particle is 
consistent with continuum dynamics. At later times, the particle is eventually adsorbed 
onto the wall and subsequently executes an intermittent stick-slip motion. We show that 
van der Walls forces are the dominant contribution to the particle adsorption phenomenon 
and that depletion forces are weak enough to allow, in the highly-wetting situation, an 
initially adsorbed particle to spontaneously desorb. 

PACS numbers: 47.15. Gf,47.11.+j,47.15.Rq,68.08.-p 
Keywords: nanochannel, molecular dynamics, suspension 



I. INTRODUCTION 

The flow of particles of different shapes suspended in Newtonian fluids at low 
Reynolds numbers and transported through different geometries has already been 
studied for over a century, and the use of a continuum description through the 
Stokes equations to model such phenomena has proved extremely successful.— There 
exist, however, situations in which the flow involves characteristic length scales far 
beyond the range of applicability of the continuum equations of fluid mechanics, and 
a microscopic, molecular detail of the fluid behavior and that of fluid-solid interfaces, 
may provide crucial understanding.— A paradigmatic example of this situation refers 
to the recurring question regarding the appropriate boundary conditions at the fluid- 
solid and fluid-fluid interfaces, a long standing problem that arises in the description 
of two immiscible fluids moving along a solid-surface.— "^ In recent years, the sustained 
development of new microfluidic devices and their potential applications in a variety 
of fields, particularly those associated with the so called lab-on-a-chip^ devices or /i- 
tas,-^ have led to a renewed interest in the problem of fluid flow under microscopic 
conflnement.— i^iiSiii Moreover, novel fabrication techniques are reaching resolutions 
in the nanometer scal o^^i^^ and nanofluidic devices are beginning to emerge as pow- 
erful tools in biotechnology and related areas.— One of the fundamental advantages 
that renders miniaturization so attractive in biotechnology is its potential ability to 
detect, transport and interrogate single particles at the molecular level. Of excep- 
tional promise are those devices that take advantage of fluid flow and suspension 
flow behavior only available at the molecular scale. It is clear then that, in order 
to realize this promise and exploit the potential beneflts of such nanodevices, it is 
crucially important to understand both the transport of nanometer size particles 
through nanochannels as well as its distinctive physics under nano- confinement. In 
this context, molecular simulations provide a powerful tool that have been success- 



fully applied to the study of the effects of nano confinement in single fluids, and 
have thereby illustrated a variety of novel and interesting dynamic behavior which 
drastically differs from the hydrodynamics of bulk fluids.-^^*^^*^ 

In a previous letter, we investigated by means of Molecular Dynamics (MD) sim- 
ulations the low Reynolds number motion of a solid sphere in a fluid-filled capillary 
tube when the dimensions of both the particle and the channel approach the molec- 
ular scale. ^'^ Specifically, we examined the behavior of the system as the wetting 
properties of the fluid are varied from perfect to partial wetting, and identified a 
novel adsorption phenomenon in which, for poorly wetting fluids, a solid sphere, ini- 
tially moving along the center of the nanochannel, is adsorbed onto the tube wall 
while displacing all the fluid molecules from the particle-wall gap. Here we present 
a detailed study of such adsorption phenomena, with particular emphasis on: i) The 
nature of the transition from adsorption to no-adsorption as the wetting improves, 
and its dependence on the shape of the particles; ii) The role and relative importance 
of the depletion and dispersion forces on the adsorption phenomena; iii) The behav- 
ior of elongated particles, which exhibit orientational selectivity upon adsorption; 
iv) The stochastic behavior of the adsorption phenomena, particularly the times at 
which adsorption takes place, the reversibility of the adsorption process, and the 
correlation between the stick-slip motion displayed after adsorption and the velocity 
fluctuations of the particle. 

II. MOLECULAR DYNAMICS SIMULATIONS 

The numerical calculations are based on standard MD techniques applied to the 
system shown schematically in Fig. specifically, a colloidal particle moving under 
the action of an external force through a periodic cylindrical channel. The particle, 
the wall of the channel, and the fluid are all treated as atomic systems interacting 
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FIG. 1: Schematic and explicit views of a colloidal spherical particle moving in a nanochan- 
nel. 



via Lennard- Jones (LJ) potentials, 
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Here r is the interatomic separation, a is roughly the atomic diameter and will be 
used as a length scale, e is the depth of the potential well and will be used as an energy 
scale, and A is an adjustable parameter which controls the attraction between the 
various atomic species. Although we consider generic non-polar materials, typically 
cr is a few Angstroms, e = O(10^A;b), where ks is Boltzmann's constant, and the 
intrinsic atomic time scale r = y^ma'^/e is a few picoseconds.— In the remainder 
of the paper, all quantities are expressed in LJ reduced units using a, e and r as 
characteristic length, energy and time scales, respectively. The potential is truncated 
at Tc = 2.5 and shifted by a linear term so that the force vanishes smoothly at 
the cutoff. Newton's equations are integrated with a fifth-order predictor-corrector 



scheme, with a time step 5t = O.OOSr. The temperature of the system is T = 1.0, 
which is maintained in the fluid by a Nose- Hoover thermostat with relaxation time^S 
Q = 10, and in the solid by velocity rescaling. 

The confining wall is composed of atoms of mass m^ = 100, tethered to fixed 
lattice sites by a stiff linear spring with constant k^ = 100. The wall sites are 
obtained from a cylindrical section of an fee lattice of length L^ and inner and outer 
radii R and Ro = R + Si/2, respectively. Here i = 1.71 is the fee lattice constant, 
corresponding to a wall number density p^ = 0.8. We refer to this channel as 
an ordered cylinder, and also consider a disordered version, where the equilibrium 
position of each atoms is perturbed by a random displacement in each direction, 
uniformly distributed in the range —£/6 < ^ < i/Q. 

The spherical nanoparticle is constructed in an analogous fashion, from lattice 
sites lying within a spherical volume of radius a = 3i = 5.13. Again we consider 
both ordered and disordered spheres, where in the latter case the atomic positions are 
perturbed randomly as is done for the wall, with — £/5 < ^ < i/5 . The atoms in the 
particle are however fixed at the lattice sites, allowing the motion of the nanoparticle 
to be computed by rigid body dynamics. The total force and torque on the particle 
are found by summing the individual atomic forces acting on all of its atoms, and 
Newton's and Euler's equations are integrated to give the particle's motion.— We 
assume here that the particle and channel wall are made of the same material and 
that their atoms interact via the LJ potential with A = 1. 

The fluid atoms are initialized on fee lattice sites with the same average density 
as the wall, interacting among themselves with the A = 1 standard LJ potential. 
The volume-average number density of the fluid is the same as that of the tube 
wall, pav = 0.8, which, at T = 1.0, corresponds to the fluid phase of the Lennard- 
Jones bulk system, slightly above the critical temperature,~ Tc = 0.937. The fluid is 
allowed to "melt" off the original lattice sites for a few r using velocity rescaling, and 
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then is equilibrated using the Nose- Hoover thermostat. A typical system has length 
L^ = 20i = 34.20, inner radius R = 12i = 20.52, outer radius Ro = 13.5i = 23.09. 
and contains 35860 fluid atoms, 459 atoms inside the nanoparticle, and 9920 solid 
molecules in the tube wall. FigureQshows a snapshot of all atoms after equilibration. 
A key feature in these simulations is the interaction between the fluid and the solid 
atoms, specified by the parameter A in Eq. (P), which provides a simple means of 
controlling the wetting properties of the fluid-solid system. Recent simulations using 
similar parameters showed that varying A from 1.0 to 0.5 leads to a transition from 
wetting to nonwetting behavior for a liquid drop on a solid substrate.-^^*^ Moreover, 
a common approximation for the contact angle in LJ systems, based on Young's 
law and a Laplace estimate for the surface energies , ^^^^^ is cos(^) = — 1 + 2A. In 
our simulations, we vary A from to 1, corresponding to estimated contact angles 
ranging from vr to 0. A second interpretation of A arises from MD simulations of 
dilute gas flow in the Knudsen regime, where^Z. A = 1 corresponds to an ideal thermal 
wall in which a molecule, colliding with the solid, undergoes a diffusive collision and 
emerges from the wall with an uncorrelated random velocity, distributed according 
to a Maxwell distribution at the temperature of the wall. The opposite limit, A = 0, 
corresponds to specular reflection where the normal velocity is reversed while the 
tangential component is preserved. 

III. FLOW OF A SINGLE FLUID THROUGH A CYLINDRICAL 
NANOCHANNEL 

The confinement of fluids at nanometer scales has a major impact on their rhe- 
ological behavior, mainly due to the increasing influence of the interfaces as the 
surface-to-volume ratio is increased. We shall show, for example, that a variation in 
the solid-fluid interaction, which induces an order a displacement in the location of 



the first peak of tlie fluid density profile adjacent to the wall (the position of the first 
fluid layer that is adsorbed onto the wall, see Fig. |2)) leads to noticeable changes in 
the average density at the center of the tube, in contrast to macroscopic flows. In 
turn, these small but noticeable changes in density induce substantial changes in the 
fluid viscosity. Thus, it is important to characterize first the behavior of single fluid 
flows in nanochannels under different wetting conditions, to provide some necessary 
background for the nanoparticle simulations. We begin, therefore, by discussing how 
the flow characteristics of the pure fluid vary with channel radius and wetting be- 
havior. We consider Poiseuille flow, in the ordered-wall case, which is generated by 
applying to each fluid atom a body force / in the axial direction. After an equili- 
bration period of 5r, this force is linearly ramped up from zero over 50r, and then 
held constant at / = 0.01 for 500r while the time-averaged density and velocity 
profiles are measured. We average over five realizations each for three values of tube 
(inner) radius, R = (4, 6, 12)£. The qualitative features of the results are the same 
for each radius, although the differences in detail merit discussion, so we show the 
profiles for the intermediate case in Fig. |21 and some numerical results in Figs. |20 
and Table m For R = 12i and A = 0.1, the average fiuid velocity is Vr^ ~ 0.3 and 
thus the Reynolds number is Re ~ 2, which is therefore the largest Reynolds number 
in our simulations. 

The density profiles show a wall-induced layering similar to that observed in liquids 
confined between planar surfaces, in both experiments^^ and MD simulations ,^^'^'^^ and 
also predicted theoretically for the case of hard-sphere systems.-^ It is important to 
note that we have also performed equilibrium simulations in the absence of fiow and 
found that the density profile agreed exactly, within the available sensitivity of the 
numerical simulations, with the non-equilibrium profiles obtained in the presence of 
fiow. It is impressive that a system driven so far from equilibrium by extremely 
large shear rates, ~ 10^° s~^, has the same statistical structure at the microscopic 
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FIG. 2: Pure fluid profiles for channel radius R = 6 £ = 10.26; open and solid circles cor- 
respond to A = 1.0 and A = 0.1, respectively. The data points are obtained as histograms 
by dividing the interior of the tube into cylindrical shells of thickness 0.05 i = 0.0855. (a) 
Density vs. r; the curves are drawn through the data points to guide the eye. (b) Axial 
velocity vs. r; the curves are best fits to a parabolic profile. 

level (singlet distribution function) as the system in thermodynamic equilibrium.— 
In most cases the density reaches a constant well-defined value at the center of the 
tube, although its magnitude differs from the volume-average density pav = 0.8, 
a consequence of the attraction or repulsion of fluid atoms near the wall. But, 
when A = 1.0, in all cases the density at the center remains essentially equal to 0.8, 
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reflecting the fact tliat tlie fluid and solid have the same density and that their atoms 
have the same interactions. On the other hand, since, even for A = 1.0, there exists 
still a layer of fluid adsorbed at the wall which slightly depletes atoms from the bulk, 
the density at the center is reduced, as seen in Table HI The number of noticeable 
layers at A = 1.0 is 5-6. In contrast, in simulations with poor wetting, A = 0.1, the 
fluid avoids the wall and the first liquid layer is displaced towards the tube center, 
resulting in an obvious increase in the fluid density at the center. A larger number of 
layers is observed in this case, 6-7, probably due to the higher packing as a result of 
the same wall-repulsion effect, and, for the smallest systems, the density oscillations 
persist all the way to the center of the channel. 

The velocity proflles in Fig. El (and in the other cases not displayed) are also 
similar to those observed in LJ fluids conflned between planar surfaces.— >2& Except 
for the region near the wall, the data are accurately fltted by a parabolic proflle 
Vxif) = ci — C2r^, from which we can obtain an effective viscosity /i = pf /Ac2 (Note 
that this way of computing the effective viscosity is different from the more standard 
calculation in which the stress tensor in computed in a simple shear flow^). In 
addition, a slip length can be deflned either as the distance from the wall where the 
fitted velocity extrapolates to zero, ^i = (ci/c2)^/^ — -R, or through a Navier condition 
relating the fluid velocity Vs and the shear stress ag at the surface, pVs = S,2(^si which 
gives ^2 = [vx/{,—dvx/dr)]\R = (ci — C2-R^)/(2c2-R). The two estimates of the slip 
length agree, as seen in Table U and in Fig. |3] As was said earlier, in the A = 1 
wetting case, there is little variation in the fluid density at the center of the tube, 
hence only small variations in the fluid viscosity are observed, with mean values 
ranging from /x = 1.83 for R = A i to fi = 2.05 for R = 12 i. On the other hand, 
larger variations in the central density for the non- wetting A = OA situation lead 
to even larger variations in the measured viscosity, which now varies from /i = 3.3 
for R = 4 i to fi = 2.63 ± 0.01 for R = 12 i. However, in all cases, the viscosity 
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FIG. 3: Shear viscosity for the Lennard- Jones fluid as a function of the central density. 
The solid circles, triangles and squares are the simulation results listed in Tabledfor R =4, 
6 and 12£, respectively. The open circles are published numerical results for the shear 
viscosity of LJ bulk systems^ and the solid line is an empirical fit to the experimental 
shear-viscosity data for Argon?^ //(p) = /^o + 0.0324(exp(5.18/9) — 1), where ^o is the 
dilute-gas value limit of the viscosity a1^ T = 1. 

variation correlates closely with the density variation, and agrees with previously 
reported experimental and numerical results for the shear viscosity of liquid Argon, 
as shown in Fig. El The slip length remains constant and negative for all tube sizes 
in the wetting simulations, ^ ~ (t/2, due to the locking of the first layer next to the 
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FIG. 4: Slip length for the Poiseuille flow in a narrow tube as a function of the wetting 
parameter A. The solid (open) squares, circles and triangles correspond to simulation 
results for the slip length ^i (^2) in tubes of radius R =4, 6 and 12 £, respectively. 

solid surface, and again is in complete agreement with previously reported results in 
planar geometries.— The non- wetting case A = 0.1 shows larger slip at the fluid-solid 
interface, as can be seen in the velocity proflles, since the fluid is no longer locked to 
the wall,— but the slip length is still 0{a), also in agreement with previous results 
obtained in planar Couette and Poiseuille flows with similar fluid densities.-^ 
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TABLE I: Numerical results for pure fluid flow. Average central density, Poiseuille fitting 
results, viscosity and slip lengths. 
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IV. MOTION OF A SPHERE THROUGH A CYLINDRICAL 
NANOCHANNEL: SHORT TIMES 

We now consider the motion of a spherical solid particle forced through a cylindri- 
cal tube of nanometer dimensions. As discussed in our previous work,^- the motion 
of the suspended Brownian particle, initially located in the middle of the tube, can 
be divided into two different regimes. The first occurs at early times, for which the 
particle has not deviated substantially from its centerline motion and, thus, does not 
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directly interact with the tube wall. This initial regime was shown to be amenable 
of a continuum description. The second regime correspond to the long time behav- 
ior of the particle, after the diffusive motion has allowed it to sample the accessible 
cross section. In this section we shall examine the behavior at early times and focus 
on the particle's mobility. The ordered spherical particle of radius a = 3 £ = 5.13 
is initially located at the center of the tube and is allowed to equilibrate with the 
fluid for lOr, whereupon a body force is ramped up from zero to / = 0.1 over a 40r 
interval. We performed simulations for tube radii ranging from R = 4 i to R = 14 i, 
corresponding to relative sizes 4/3 < R/a < 14/3, and for a length of the (ordered) 
tube L^ = 30 i. 

Periodic boundary conditions are used in the axial direction, implying that the 
simulations in fact correspond to the motion of the "real" sphere and an infinite 
string of periodic replicas. Since hydrodynamic interactions between suspended solid 
particles are usually long-ranged, one might be concerned about the interactions be- 
tween the periodic images. Fortunately, however, confining geometries effectively 
screen such hydrodynamic interactions and lead to a much faster spatial decay of 
the disturbances created by a single particle, typically exponentially with a charac- 
teristic length fixed by the tube radius^ R. In particular, the interaction between 
spheres, equally spaced along the centerline of a tube, has been shown to be small 
when the separation between spheres is larger than one tube diameter,'^ and re- 
cent experiments involving colloidal particles in a narrow channel have demonstrated 
that hydrodynamic interactions are indeed strongly inhibited by the one-dimensional 
confinement.-^ Given that, even in the worst scenario, we have L^/R > 2, we can 
then safely neglect any hydrodynamic interactions with the periodic copies of the 
particle. 

The radial motion at short times is indicated in Fig. for A = 1.0 and for three 
different values of the tube radius. The motion is roughly random and diffusive, but 
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FIG. 5: Radial position of the spherical particle, a = 3 £, as a function of time for tube 
radii R = A i (open circles), R = 6 £ (open triangles) and R = 14 i (open squares). Time 
t = corresponds to end of the force ramp. In all cases there is perfect wetting, ^4 = 1.0. 

its range in the radial direction is constrained by the presence of the tube, decreasing 
as the tube radius approaches that of the particle. To examine the axial velocity, 
we follow the evolution of the particle for 300r, measuring the velocity every 2r 
and averaging the last 100 points (corresponding to a time average over about half 
a nanosecond). In Fig. IHl we present the average axial velocity of the sphere, Ux, 
as a function of the relative size of the tube R/a. Note the large error bars in the 
figure which reflect the substantial fluctuations in the velocity of the particle due 
to the thermal motion. The axial velocity is normalized by the Stokes velocity, the 
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FIG. 6: Particle average axial velocity, normalized by the Stokes velocity Ug (Hav = 1-94; 
Na = 438), as a function of the relative size of the tube radius: A = 1.0 (closed circles), 
A = 0.5 (open triangles) and A = 0.25 (closed squares). The solid curve is the theoretical 
prediction of Ref. |4l|, as discussed in the text. 



velocity which the particle would have had in an unbounded fluid under the action 
of the same external body force, Us = F/lGnfia), where F = Naf is the total force 
acting on the particle and Na the total number of atoms in it. Using the measured 
average viscosity as described in Section IIIIl /!„„ = 1.94, and the average number 
of atoms in a spherical particle, Na = 438, we obtain Ug = 0.23. (A similar value. 
Us = 0.22, is obtained using the experimental value of the viscosity at the volume- 
average density and the average force on a sphere with the same density). There 
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FIG. 7: Particle average axial velocity as a function of the wetting properties of the system 
(The tube radius is i? = 6 ^ and the particle radius is a = 3 ^). (a) The velocity is given 
in MD units, (b) The velocity is normalized by a Stokes velocity U* that accounts for the 
variations in the viscosity, U* = F/(6irfi*a), with /x* being the viscosity values listed in 
Table n 

is some ambiguity in the Stokes velocity arising from an ambiguity in the particle 
radius in that, if instead of taking the radius of the suspended particle as a = 5.13, we 
were to define the effective hydrodynamic radius as that of the volume inaccessible 
to the fluid molecules,— aes ~ a + 1, the predicted Stokes velocity would decrease 
by about 20%. This uncertainty in the definition of the particle radius is intrinsic 
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to its molecular nature and the small length scales involved in the problem. Also 
note that the Stokes velocity, Us = 0.23, was computed using the average viscosity 
obtained for A = 1.0 and it does not account for the variations in the viscosity that 
occur due to changes in the wetting properties (see Fig. ^. As shown in Fig. Efa), 
there is little statistically significant variation of the velocities with changes in the 
wetting parameter A, particularly below about A = 0.5. On the other hand, when 
the velocity is normalized by a Stokes velocity that accounts for the variations in 
the viscosity, U* = F/iQ-r^jj^a), with /i* being the viscosity values listed in Table HJ 
a completely different picture emerges, as can be observed in Fig. IZfb), showing a 
systematic increase in the velocity with decreasing A. We can then conclude that 
the apparent lack of dependence of the average velocity on the wetting properties, 
observed in Fig.[7|^a), is in fact due to the competition between an increase in the fluid 
viscosity and a reduction in the drag coefficient of the particles as A is decreased. 

Next, we wish to test the degree to which these results can be described by con- 
tinuum hydrodynamics. Previous studies have shown good agreement in a number 
of situations,— including the Stokes drag on a sphere in an unbounded fluid,-^ even 
where the particle size is comparable to that of the liquid molecules, and the motion 
of a spherical particle through a cylindrical pore^^ where the tube and particle have 
molecular structure and dimensions similar to those in the present work. (Unfor- 
tunately, the particle's axial mobility was not considered in R,efl44l) The particle's 
radial Brownian motion and molecular-scale surface roughness are difficult to address 
in standard continuum calculations, so we compare our results with the solution of 
the Stokes equations for a non-Brownian smooth sphere moving along the center of 
a smooth tube. Note that although the particle Reynolds number is Re ~ 0.5, com- 
puted using the Stokes velocity of the particle in an unbounded fluid, inertia effects 
will not be accounted for in the continuum calculation. The continuum description 
we used also neglects density variations across the tube, and the associated viscos- 
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ity variation discussed in Section IIIIl It is also assumed that the no-shp boundary 
condition apphes at both sohd-fluid boundaries and thus we expect this continuum 
description to more accurately describe the case A = 1.0 for which slip effects were 
shown earlier to be negligible in the pure fluid problem in accordance with previous 
studies.-^"^ 

The velocity of the particle is determined only in part by the force applied to 
it, since one must also specify some additional information about the motion of the 
fluid, such as the pressure drop along the channel, AP, or the mean fluid velocity 
along the tube, Vx- Actually, it is convenient to work with the excess pressure drop 
due to the presence of the sphere, AP, = AP — APq, where APq = SfiL^Vx/R"^ is 
the Poiseuille pressure drop for the pure fluid. The hnearity of the Stokes equations 
and the no-slip boundary condition imply that 



F 

fvrp2)AR 



fia 



TZfu T^fv 
TZpu TZpv 



Ux 

Vx 



(2) 



where the resistance matrix elements TZfj,i, are functions of R/a alone when the par- 
ticle moves along the center of the tube, and were given explicitly by Bungay and 
Brenner .i^i*^ 

In the case at hand, we have periodic boundary conditions in the x direction and 
the total pressure drop is zero. Therefore {nR'^)APs = —{7iR'^)APq = — (/ia)7^o'^x, 
where T^o = SnLx/a. Equation (0) can be solved for F and Vx in terms of Ux, and 
in particular F = {fia)TZUx where 

TZfvT^pu 



n = Ufu - 



(3) 



7?.o + T^pv 

In Fig. ini we compare the velocity, as predicted by the above expression, to the 
numerical results, and see that the agreement is satisfactory, given the substantial 
error bars (due to thermal fluctuations). Surprisingly, the numerical results are 
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not sensitive to variations in the wetting properties of the system, with a lack of 
apparent trend in the average velocity as A is decreased from A = 1.0 to A = 0.25, 
and the continuum approximations seem to describe all possible wetting conditions. 
However, note that, by normalizing the particle velocity by a single Stokes velocity, 
which was computed using the average viscosity fiav obtained for A = 1.0, we have 
not eliminated the dependence on /x. In particular, its dependence on the wetting 
parameter A is not accounted for in Ug. Thus, we can conclude again that, owing 
to the competition between an increase in the viscosity as A is decreased, and the 
concurrent reduction of the drag on the particles due to a larger slip-effect, the results 
are strikingly insensitive to variations in the wetting properties of the fluid. 

One can show, from Eq. (jH} and the explicit forms of the resistance matrix el- 
ements given in Ref. |4l|, that the continuum mobility vanishes as (1 — R/aY^'^ as 
a ^ R, but, although this behavior is consistent with the simulations, it does not 
appear to capture the trend. Indeed, when the gap between the particle and the wall 
approaches a molecular diameter, one might expect a continuum result to fail, but 
we do not have precise enough statistics to make a definitive statement. 

V. MOTION OF A SPHERE THROUGH A CYLINDRICAL NANOCHAN- 
NEL: LONG TIMES 

The results presented thus far for the mobility of the particle correspond to rela- 
tively short-time motions, essentially at the center of the channel, where any direct 
interactions with the wall can be ignored. After a sufficiently long time, however, 
the thermal motion of the particle should cause it to sample the whole cross-section 
so that direct interactions with the tube wall are unavoidable. An estimate of the 
(diffusive) time required to reach the wall is r^, r^ {R — ay/2D, where the diffusion 
coefficient in turn can be estimated from its bulk value given by the Stokes-Einstein 
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FIG. 8: Radial position of the sphere vs. time for a single realization at four values of A: 
0.1, 0.4, 0.7 and 1.0. The sudden jump to a radial postion ~ 5 at the three lower values of 
A corresponds to adsorption at the tube wall. (The tube radius is i? = 6 ^ and the particle 
radius is a = 3 i. The dashed line corresponds to a radial postion equal to R — a, the 
nominal maximum radius.) 

formula, Z^buik = ksT/GiTfia. We find td ~ 2500r, roughly 10 times the elapsed time 
in the simulations described earlier. Therefore, in order to explore wall interactions, 
we have conducted rather longer simulations, O(lO^r), in a channel of radius R = 2a 
for a range of values of A, for a disordered particle and wall. 

Figure IHl illustrates the effects of the wall interaction, showing the radial position 
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of the particle as a function of time for four independent realizations in which the 
fluid changes from non- wetting to wetting, corresponding to A = 0.1, 0.4, 0.7 and 
1.0. The other parameters of the simulation are held fixed: R = Gi; a = 3i, Lx = 30i 
and / = 0.1. It is seen that, except for the completely wetting case, the particle 
eventually reaches a critical distance from the wall, and then rapidly moves to the 
wall and becomes adsorbed, with the adsorption time decreasing with A. Note 
that, because the sphere and the tube are disordered and the LJ interaction is soft, 
a particle can reach a radial position with its center located beyond the nominal 
maximum value R — a = 5.13. The snapshots of the final configuration in the 
simulations in Fig. IHl show that the adsorbed particles are in direct contact with the 
wall, having expelled the fluid from the gap between them, except for the completely 
wetting case where an intervening fluid layer persists. To be sure, in any continuum 
description to this problem, a lubrication force resisting the particle motion would 
be generated which would diverge as the gap between particle and wall vanishes.-^ As 
in earlier MD calculations, however, this divergence cuts off at molecular distances, 
and in fact the resistance force vanishes when all of the fluid drains out of the gap.-^ 
We can obtain some insight into the particle adsorption by considering the problem 
in terms of van der Waals interactions between solids, represented by the attraction 
term of the LJ potential — 4e/r^. In a continuum description, the attraction energy 
between a solid wall and a solid sphere of radius R, separated by an empty narrow 
gap of thickness h '^ R, is given by^i^^S*^ 

>4.W = -^4(l-^)"', (4) 

if the interaction is non-retarded and additive, where Ah = Aen'^psPw is the Hamaker 
constant of the system.— "^ (Note that in MD simulations both assumptions are 
implicit in using a constant interatomic potential.) The energy in this "Derjaguin 
approximation" exhibits limitations similar to those of the continuum description 

22 



when applied at molecular scales in the case of the lubrication forces, in that the 
expression for Vgw given above diverges at contact. In fact, since the short range 
repulsion present in the L J potential would prevent molecules from approaching closer 
than a separation distance of order a, we could use Vsw{<^) as a reasonable estimate 
of the energy at contact---^ (Note however that the snapshots in Fig. El show that the 
actual wall-sphere configuration corresponds to a more complex inter-meshing of the 
two solids, which would be rather difficult to capture in any continuum description.) 

We can correct Eq. (JD) for the effect of the suspending fluid, by taking advantage 
of the additivity of the interatomic potential in MD simulations. Specifically, to the 
expression for Vsw{h) given above, we can add a term accounting for the interaction 
of the tube wall with the fiuid completely filling the interior of the tube, from which 
we then subtract the interaction of the tube wall with a sphere of fluid, located at the 
same position as the solid particle. The former interaction potential is a constant, 
independent of the position of the solid sphere and can therefore be ignored, while 
the latter interaction term is of the same form as Eq. (^ but multiplied by the 
wetting parameter A. Thus, the expression for the potential remains the same as in 
Eq. Q, but with an effective Hamaker constant given by A'f^ = 4e(l — A)tt'^PsPw. 

The corresponding attractive force between the sphere and the tube wall is given 
by the derivative of Eq. Q with respect to the separation distance h, 

which is proportional to {1 — A). In Fig. ITTH we show that the adsorption transition 
strongly depends on A, in that the sudden jump into contact observed for A = 
0.1, becomes less abrupt as the fluid becomes more wetting, and 1 — A decreases, 
which is qualitatively consistent with the previous result that the attractive force is 
proportional to 1 — A. 

In addition to these energetic considerations for particle adsorption, another pos- 

23 



sible inechanism, leading to attraction between solid bodies in a fluid, is a depletion 
fQj,pQ48,49,50,5ij52 j-gg^lting from the reduction of the excluded volume and the conse- 
quent increase in the entropy when solid particles are in close proximity. However, the 
results given previously show that changes in the interaction energy have a dramatic 
effect on adsorption, and suggest that the dominant considerations are energetic as 
further indicated from the following "inverse" simulations. Specifically, we consider 
the evolution of an initially adsorbed colloidal particle suspended in a non-wetting 
fluid, A = 0.1, when the fluid is instantaneously rendered completely wetting by 
changing A to 1.0. After this change, the intermolecular forces between all species 
are identical, so one might expect that these forces are approximately isotropic, and 
that only depletion forces would be available to preserve the adsorption of the par- 
ticle. In Fig. El we present the evolution of the radial position of the sphere in two 
independent inverse simulations. In both cases, the wetting parameter is initially set 
at A = 0.1 which, in time, leads to the sphere being adsorbed onto the tube wall; 
then, after allowing the sphere to equilibrate in an adsorbed position for t = lOOOr, 
the wetting parameter is switched to A = 1.0. As is seen in Fig. ^2 the particles 
desorb spontaneously with no measurable trapping time, thereby leading us to con- 
clude that depletion forces are not responsible for the long (perhaps irreversible) 
adsorption times observed in the poor wetting cases. This calculation is of course 
somewhat indirect in that we do not estimate directly any entropic effects, but in a 
forthcoming paper^^ we calculate the Helmholtz free energy changes associated with 
particle motion in nanotubes, which we hope will shed some light on this question. 
Finally, we discuss the behavior of particles a/ter adsorption. In Fig.^Jwe present 
the radial and axial positions of a disordered sphere as it moves through a disordered 
tube, for the case A = 0.2. It is clear that, even after adsorption has taken place, the 
particle continues to move along the tube in the form of an intermittent stick-slip 
motion while remaining in near contact with the tube wall. This form of a stick-slip 
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motion has no counterpart in any simple continuum description, because, even if it 
were possible for the sphere to overcome the resistance of the lubrication layer of 
the fluid and almost touch the cylindrical wall, say at a separation h <^ a, the force 
required to move it along the surface would diverge as^^ a/h, which far exceeds our 
applied force. Other types of interaction, such as double-layer effects, could lead to 
an equilibrium radial position of the sphere at some distance away from the wall,— 
but this also would not account for the stick-slip behavior observed here. 

Stick-slip motion was found to occur for all types of ordered and disordered spheres 
and walls, and thus does not depend on any precise matching of their underlying 
molecular structures. In fact, the results shown in Fig. El correspond to the extreme 
case where the position of both the molecules of the sphere and those of the tube 
wall have been randomly modified. In our previous paper— we reported that, both 
the average radial position of the sphere and the random fluctuations around that 
average, depend on whether the particle is sliding along the wall {slip) or whether 
it is motionless on average {stick). For the probability density function of the radial 
position of the solid particle, we found similar results for different types of parti- 
cles and tube wall, but the distributions corresponding either to sliding or sticking 
particles have substantial overlap in the case of disordered solids, and the average 
distance to the wall in both cases is almost indistinguishable. On the other hand, in 
all cases there exists a strong correlation between the amplitude of the fluctuations 
in the radial position and the motion of the particle in the axial direction. This 
is illustrated in Fig. El which shows how, during a single realization in which the 
particle displays intermittent stick-slip motion, the amplitude of the fluctuations in 
the radial position of the particle is clearly correlated to its axial motion. 
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VI. MOTION OF A PROLATE SPHEROID THROUGH A CYLINDRICAL 
CHANNEL 

In practice, one may be interested in particles of more general shapes, and we now 
consider the motion of prolate spheroids through a cylindrical nanochannel. These 
spheroids were constructed in a fashion analogous to the spherical solid particles 
discussed in section |TT1 with the minor axis equal to a, the radius of the spherical 
particles considered previously, and the major axis b related to the ellipticity e via 
e = a/1 — {a^ /IP'). The simulations were performed for two different aspect ratios, 
a = b/a, a = 3 corresponding to e = 0.89 and a = 5 corresponding to e = 0.96, and 
both ordered and disordered spheroids were studied. The tube radius always equaled 
twice the minor radius of the spheroids, R = 2a, and its length was changed to 
Lx = 50£ for the case of aspect ratio a = 5, in order to avoid possible hydrodynamic 
interactions between periodic images. 

Initially the spheroids were placed at the center of the tube with their major axes 
oriented along the tube axis. Note that neither of the two spheroids would have fitted 
inside the tube if its major axis had been oriented perpendicular to the tube axis, 
with the maximum possible (polar) angle between the major axis of the spheroid 
and the tube axis being Omib/a = 3) < 7r/4 and Omip/a = 5) < 7r/8. As in the case 
of spherical particles, the surrounding fluid was first equilibrated and then a body 
force, linearly ramped up from zero to a final value / = 0.1, was applied to each atom 
in the spheroid. We followed the evolution of the spheroids as they moved through 
the tube for relatively long times, up to lO^r ~ 0.25/is in some cases. 

The qualitative behavior of the spheroids is very similar to that of the spheres 
under similar conditions. Specifically, in section we showed that the adsorption 
events are readily observable in the evolution of the radial position, as revealed 
by a sudden jump-into-contact of the particles. In Fig. El we present the time 
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evolution of the polar angle 9 for independent MD simulations corresponding to 
different wetting conditions, varying from complete wetting A = 1.0 to non- wetting 
A = 0.1. As in the case of the sphere, for poor wetting (low A) the polar angle 
suddenly locks at a certain value, a transition that corresponds to the adsorption of 
the spheroid, while, in the complete wetting situation, the spheroid does not adsorb 
even after relatively lengthy simulations. Completely analogous results are obtained 
for a spheroid having a larger aspect ratio, a = 5. In addition, and going beyond 
the observations for spherical particles, the spheroids display a certain selectivity of 
the orientation angles at which adsorption takes place, and this orientation remains 
essentially locked after adsorption. In some cases, however, a transition between two 
different orientations can be observed after the particle has been adsorbed (see Fig.HSl 
for the A = 0.4 case). In Fig. El we present the probability distribution function for 
the orientation of the spheroids after the occurrence of adsorption, obtained from 
more than 50 different realizations for spheroids with aspect ratio a = 3, and with 
the wetting parameter varying from A = 0.1 to A = 0.7. We see that three different 
polar orientations predominate, corresponding to: i) spheroids with their major axis 
parallel to the tube, 6 sin 6* ^ 0; ii) spheroids occupying one fourth of the tube cross 
section 6 sin 6* ^ -R/2; and iii) spheroids locked across the tube 6 sin 6' ^ R. Snapshots 
of the three possible orientations are presented in Fig. ^J for the case of disordered 
spheroids and a disordered tube. We can also see in this figure the complete depletion 
of fluid molecules from the narrowest gap between the tube wall and the spheroid. 

Finally, we also observe that, after adsorption, the spheroids display, in some 
cases, a stick-slip motion similar to that observed for adsorbed spheres. In Fig.lT^we 
present the time evolution of the radial position and of the orientation for an adsorbed 
solid spheroid with a = 3 and A = 0.5. Clearly, there exists a strong correlation 
between the amplitude of the fluctuations in the polar angle and the motion of the 
spheroid in the axial direction, in that, larger angular fluctuations correspond to a 
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sliding motion, while a sticking behavior corresponds to much smaller fluctuations 
in the orientation. 

VII. CRITERION FOR THE ADSORPTION TRANSITION 

In the previous sections we showed that both spherical and elongated particles 
suspended in a poorly wetting fluid are eventually adsorbed onto the tube wall. In 
the simulations, the wetting properties are controlled by the parameter A in Eq. (fT|). 
and a simple estimate of the attractive force between the particle and the wall in 
Eq. (0) likewise involved A, so it is natural to ask whether there exists a critical value 
of A which identifles an adsorption transition. To this end, we simulated more than 
200 cases of spheroids of aspect ratios 1, 3 or 5, in fluids with different values of A, 
and with either ordered or disordered atomic structures for the particle and for the 
tube wall, with multiple realizations in some cases. 

Table |n] summarizes the parameters of the various runs (all with R/a = 2). 
The simulations ran for long times, O(lO^r) or longer, and we can characterize the 
results numerically in terms of the average particle velocity, with the result shown in 
Fig. El In this flgure, the velocities are normalized by the Stokes velocity that the 
same particle would have had in an unbounded fluid with the same applied force,— 
U,{a) = F/{6nfxaK{a)), where K{1) = 1, K{3) = 1.40446, and K{5) = 1.78481. 
Relative velocities of 0.1 or below correspond to the stick-slip motion of adsorbed 
particles already described, while values near 0.3, which are similar to the velocities 
obtained at short times (c/. Fig. IHlat R/a = 2), correspond to the motion of spheres 
that are not adsorbed. 

There exists strong evidence that the transition to adsorption occurs ai A ^ 
0.8 independent of particle characteristics. Moreover, we also performed additional 
"inverse" simulations to those discussed in Section IVl We found that whereas, for 
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TABLE II: Adsorption times in MD simulations. Note that an estimate of the diffusive 
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P: Particle type, D Ordered; jZl Disordered. 
C: Tube wall type, D Ordered; jZ] Disordered. 
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Tt: Duration of the simulation in units of lOOOr. 



A > 0.8, the particles are soon desorbed from the wall, for A < 0.8 the particles 
remained adsorbed for the duration of our simulations (e.g. lO^r for A = 0.6), 
thus providing further evidence for the existence of a sharp transition. Since particle 
adsorption is a random process which requires a fluctuation to bring the particle near 
the wall, the time required can vary over a wide range (see table IVTT|) . and there is 
no guarantee that the transition to adsorption will occur during the duration of the 
simulation. Indeed, as is seen in Fig. El there is one exception to the general rule, 
for a sphere a.t A = 0.6. 

VIII. SUMMARY AND CONCLUSIONS 

The complex dynamics of suspended nanometer size particles under nano- 
confinement was investigated by means of molecular dynamics simulations. We 
studied the low Reynolds number transport of nanometer size spheroids with ra- 
dius a = 3i, and aspect ratios ranging from 1 to 5, suspended in a L J fluid of density 
p = 0.8 and confined inside a nanochannel of radius ranging between R = A£ and 
R = 12L 

We initially examined the behavior of single fluids flowing through a nanochannel 
for varying wetting conditions, and showed that the increase in surface-to-volume 
ratio, typical at the nanometer scales, changes the transport properties of the fluid 
in a significant way. Specifically, we showed that a change in the strength of the 
solid-fluid interaction leads to a small shift in the position of the first adsorbed layer 
of the fluid, < a, which in turn leads to a relatively small increase in the fluid density 
at the center the tube, Ap/p < 13%. This small change in density, however, has 
a significant effect in the corresponding fluid viscosity at the center of the tube, 
which shows an increase of up to 80%. However, this dependence of the viscosity on 
the density was shown to be in agreement with that observed in bulk fluids. The 
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measured slip lengths are also in agreement with measurements obtained in bulk, 
or in substantially larger geometries. The structure developed by the fluid close to 
the tube wall is also consistent with previous flndings as well as with equilibrium 
simulations. Finally, we showed that macroscopic continuum equations that allow 
for slip at the solid boundaries accurately describe the velocity proflle of a fully 
developed flow in a pipe, even in the region where the fluid is highly structured. 

We then considered the motion of a spherical solid particle through a nanochan- 
nel, at both long and short times compared to the diffusive time required for the 
sphere to reach the wall. We flrst showed that, at short times, the average velocity 
of the particles remains essentially unaffected by changes in the wetting parameter 
A, particularly so for A < 0.5, and concluded that this lack of dependence on A of 
the mobility of the particles is due to the competition between an increase in the vis- 
cosity of the fluid and the reduction of the drag force on the particles as A decreases 
and the degree of wetting is reduced. Then, we compared the numerical results for 
the mobility of the spheres with the prediction based on continuum calculations for 
a non-Brownian smooth sphere moving along the centerline of a structureless tube in 
the limit of vanishingly small inertia. We found reasonable agreement, especially in 
view of the simplifled continuum model that we used which does not account either 
for the large thermal fluctuations present in the system, or for the molecular struc- 
ture of both solids. Moreover, due to the competition between drag reduction and 
increased viscosity as the degree of wetting is reduced, the continuum representation 
describes fairly well the numerical results for all wetting conditions, even though the 
assumed no-slip boundary condition becomes less accurate as A decreases. At long 
times, however, the interactions of the suspended particle with the tube wall become 
important and we showed, that for poorly wetting suspending fluids, A < 0.8, the 
particles eventually are adsorbed to the tube wall while, surprisingly, displacing all 
fluid atoms from the gap, and, subsequently, exhibit stick-slip motion along the tube 
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wall. These interesting features could never have been captured by any simple con- 
tinuum description, where lubrication forces would prevent the particle from either 
moving towards the wall or along it for small separation distances {h <^ a). We 
then studied the contribution of van der Waals and depletion forces to the particle 
adsorption phenomena and concluded that the dominant contribution comes from 
the fluid-solid van der Waals interactions, as suggested by the dramatic effect that 
changes in the wetting parameter A have on the adsorption process. Simulations 
showing the spontaneous desorption of particles in the complete wetting situation, 
A = 1.0, confirmed the conclusion that depletion forces are weak enough to allow for 
the desorption of particles over extremely short times. 

The behavior of prolate spheroids of two different aspect ratios, a = 3 and a = 5, 
is qualitative similar to that of the spheres, in that adsorption of the particles takes 
place for poor wetting conditions and that a similar stick-slip motion along the tube 
is also subsequently observed in some cases. In addition, spheroidal particles exhibit 
a certain selectivity in their final orientation upon adsorption by seemingly preferring 
three different orientations, two of which appear to locally maximize the contact area 
between the particle and the tube wall while the third appears to correspond to the 
dynamic locking of the spheroids across the tube section. 

Finally, we investigated the adsorption transition as the wetting properties vary 
from non-wetting to complete wetting, and found strong evidence suggesting that 
the transition occurs a.t A ^ 0.8 independent of particle shape. A detailed study of 
the Helmholtz free energy of the system for different wetting conditions is currently 
being performed and will be the topic of a forthcoming paper.^^ 
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FIG. 9: Cross-section snapshots of the particles at the end of the simulations in Fig. |H1 
For visual clarity, we use circles of different sizes to represent the different species, with 
radii in the ratio 1.0 to 0.8 to 0.4 for the particle, wall and fluid atoms, respectively. (The 
tube radius is R = 6 £ and the particle radius is a = 3 i.) 
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FIG. 10: Radial position of the sphere as a function of time for six single realizations 
to A = 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, and 0.7. The time axes are shifted so as to be able to 
observe the jump into contact in all cases. (The tube radius is R = 6 i and the particle 
radius is a = 3 £.) 
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FIG. 11: Radial position of the sphere as a function of time for two independent realizations 
of the MD simulations. In both cases the wetting parameter is set, initially, equal to 
A = 0.1, leading to the adsorption of the particles. At t = lOOOr the wetting parameter is 
changed to ^ = 1.0, and a spontaneous desorption of the spheres can be observed. After 
changing, at t = 2000r, the wetting parameter back to ^ = 0.1 the spheres are re-adsorbed. 
(The tube radius is R = 6 i and the particle radius is a = 3 i.) 
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FIG. 12: (a) Position of the particle along the tube as a function of time (b) Radial 
position of the sphere as a function of time. The results correspond to a disordered sphere 
moving through a disordered tube for A = 0.2. (The tube radius is R = 6 i and the particle 
radius is a = 3 ^.) 
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FIG. 13: Angle between the spheroid major axis and the tube axis as a function of time 
for four realizations at ^4 = 0.1, 0.4, 0.7 and 1.0 (for each value of A one independent MD 
simulation is presented). The abrupt locking of the polar angle observed at A = 0.1, 0.4 
and 0.7 corresponds to the adsorption of the spheroid on the tube wall. The tube radius 
is R = 6i. The spheroids have minor axis a = 3£ and aspect ratio a = 3, and both the 
spheroids and the tube are ordered. 
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FIG. 14: Probability distribution function for the orientation angles of a spheroid after 
adsorption, averaged over 50 or more realizations with wetting properties ranging from 
A = 0.1 to 0.7. The left (right) hand plots are the probability distribution functions for 
the polar angle 9 (azimuthal angle (p), and the top (bottom) plots correspond to an ordered 
(disordered) spheroid and wall. The tube radius is R = Qi. The spheroids have minor axis 
a = 3i and aspect ratio q = 3. 
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FIG. 15: Three cases of spheroid adsorption in three distinct realizations for A = 0.1, with 
side views on the left and cross-sectional views on the right. The tube radius is R = 6i. 
The spheroids have minor axis a = 3£ and aspect ratio a = 3. 
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FIG. 16: (a) Position of the spheroid along the tube as a function of time, and (b) polar 
angle of the spheroid as a function of time, for a disordered spheroid moving through a 
disordered tube, for A = 0.5. The tube radius is -R = QL The spheroids have minor axis 
a = 2>l and aspect ratio a = 3. 
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FIG. 17: Long-time average particle velocity along the tube. The velocity is normalized 
by the Stokes velocity Ug that a sphere having the same radius would have had in an 
unbounded fluid, for the same applied force and no-slip boundary conditions. The time span 
of the simulations, Tt > lOOOOr, was always longer than the diffusive time required for the 
particle to reach the wall, td ~ 2500r. The average velocity of particles that were adsorbed 
during the simulation is computed for times after adsorption had occurred, whereas, for 
particles that are not adsorbed the velocity is averaged over the entire simulation. The 
dispersion in the velocities of adsorbed particles is a consequence of the stick-slip motion 
after adsorption that has already been described. 
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